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Influences of nonlinear screening on efi'ective interactions between spherical macroions in charged 
colloids are described via response theory. Nonlinear screening, in addition to generating effective 
many-body interactions, is shown to entail essential corrections to the pair potential and one-body 
volume energy. Numerical results demonstrate that nonlinear effects can substantially modify effec- 
tive interactions and thermodynamics of aqueous, low-salt suspensions of highly-charged macroions. 
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Colloidal particles, typically nanometers to microns in 
size, commonly acquire electric charge in solution [Q. 
Practical examples are latex microspheres in paints, ionic 
micelles in detergents, and platelets in clay suspensions. 
Beyond traditional applications, e.g., in the chemical, 
pharmaceutical, and mining industries, the rich thermal 
and optical properties of charged colloidal crystals have 
inspired recent proposals for novel technologies, such as 
fast optical switches and photonic band-gap materials . 

Among the various interparticle interactions at play in 
charged colloids, electrostatic interactions strongly affect 
thermodynamic stability and materials properties 
The classic theory of Derjaguin, Landau, Verwey, and 
Overbeek (DLVO) describes the bare Coulomb inter- 
actions between pairs of macroions as screened by sur- 
rounding microions (counterions and salt ions) in solu- 
tion. In recent years, the general validity of the DLVO 
screened-Coulomb pair potential has been questioned in 
light of reported observations of phase separation ||] and 
related anomalies in strongly deionized suspensions. 

The DLVO pair potential was originally derived by lin- 
earizing the mean- field Poisson-Boltzmann (PB) equa- 
tion for the electrostatic potential around a macroion [Q . 
Density-functional and linear response ^ theories 
offer alternative routes to effective interactions, but also 
yield, as an important by-product, a one-body (volume) 
energy that contributes to the total free energy. Com- 
mon to all approaches is the simplifying assumption that 
the microion densities respond linearly to the macroion 
charge density. While linearization is often a reasonable 
approximation, its justification is uncertain in cases of 
high macroion charge and low salt concentration (ionic 
strength), precisely those conditions under which anoma- 
lous phase behavior is observed. 

The purpose of this paper is to define the range of 
validity of the linearization approximation by extending 
the DLVO theory beyond linear screening. Applying re- 
sponse theory to the primitive model of charged colloids, 
leading-order nonlinear corrections to the effective pair 
potential and volume energy are obtained and shown 
to qualitatively affect thermodynamics of concentrated, 
deionized suspensions of highly charged macroions. 

The model system comprises Nm spherical macroions. 



of hard-core diameter a and uniformly-distributed sur- 
face charge —Ze (e being the elementary charge), and Nc 
point counterions of charge ze dispersed in a symmetric 
electrolyte in a volume V at temperature T. The elec- 
trolyte consists of Ns point salt ions of charge ze and Ns 
of charge — ze in a uniform solvent, characterized entirely 
by a dielectric constant e. The microions thus number 
^ Nc + Ns positive and iV_ = Ng negative, for a to- 
tal of iV^ — Nc+2Ns. Global charge neutrality constrains 
macroion and counterion numbers via ZN^ = zNc. 

Initially ignoring salt, the Hamiltonian decomposes 
into three terms: H = Hm + He + H„ic. Here Hm is 
the bare macroion Hamiltonian involving Coulomb in- 
teractions, Vmm{r) = Z'^e'^/er, between macroion pairs 
at center-center separation r > a; He is the counte- 
rion Hamiltonian involving Coulomb pair interactions, 
Vccif) — z^e^/er; and H^c is a sum over macroion- 
counterion pair interactions of the form 
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As the interaction for r < a/2 is arbitrary, a is chosen 
to minimize counterion penetration inside the core |7| . 

The mixture of macroions and counterions is for- 
mally reduced to an equivalent one-component system of 
"pseudo-macroions" by performing a restricted trace over 
counterion coordinates, keeping the macroions fixed. De- 
noting counterion and macroion traces by ( )^ and ( 
respectively, the canonical partition function is 
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where i?e// — Hm + Fc is the effective one-component 
Hamiltonian, (3 = l/fcsT, and 



(^e^p[-l3{Hc + Hmc)]) 



(3) 



is the free energy of a nonuniform gas of counterions in 
the presence of the macroions. Adding to He (and sub- 
tracting from Hmc) the energy of a uniform compensating 
negative background, Eb, converts He to the Hamiltonian 
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of a classical one-component plasma (OCP) of counteri- 
ons. In practice, the OCP is weakly correlated, with cou- 
pling parameter T = Xb/o-c *C Ij where Xg — /3e^/e is 
the Bjerrum length and = (3/47^^)^/^. The OCP has 
mean density Uc = Nc/V{l — 77), since the counterions 
are excluded (with the background) from the macroion 
core volume fraction, 77 = ^{Njn/V)(j'^ . 

The counterion free energy is approximated by regard- 
ing the macroions as an "external" potential for the OCP 
and invoking perturbation theory p, nO|: 
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where Fqcp = -ksTln (exp[-/3(i7c + Eb)])^ is the OCP 
free energy, H'^^ = Hmc — Eh, and the A- integral charges 
the macroions. In terms of Fourier components. 



(H. 



rnclx = — ^W„c(fc)/5™(k) (/7c(-k));^, (5) 



where Vmc{k) is the Fourier transform of Eq. (|l|) and 
/5m (k) and /5c (k) are Fourier components of the macroion 
and counterion number densities. Now, to second order 
in the macroion density, we may write (for k ^ 0) 

/6,(k)= X (iHfc)«mc(fcK.(k) + :^5]x^'^(k',k-k') 

k' 

X ■e„c(fc')''^mc(|k-k'|)/0™(k')/5™(k-k') (6) 

where x'"^'' and x''^^ £^re, respectively, the linear and first 
nonlinear response functions of the OCP 

To specify the response functions, we adopt the ran- 
dom phase approximation (RPA) , which equates the two- 
particle direct correlation function of the OCP to its ex- 
act asymptotic limit: c^^^(r) = —f3vccir) Within 
the RPA, x*-^^ and x*-^' have analytic forms 



where k = \/ 'inncz'^ Xb , and 
X(2) = -^'"-/^ (8) 

^ (l+K2/fc2)(l + /,-2//c'2)(l + /t2/|k + k'|2) ^ ' 

With the response functions specified, the counterion 
density can be explicitly determined from Eq. (^ . Pene- 
tration inside the macroion cores is reduced, to at most a 
few per cent, by choosing a = na/ {2+Ka) in Eq. (|l|). Salt 
is incorporated via additional response functions for the 
multiple microion species. In the process, k is modified 
by replacing the counterion density, rij,, by the total mi- 



croion density, = 



where n± — A^±/y(l — rj) 



are the densities of positive and negative microions. 

Combining Eqs. (U)-(|8|), the effective Hamiltonian can 
be written as a sum of one-, two-, and three-body terms. 
The one-body term, or volume energy, is the sum of 
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which is the linear response volume energy and 
AS = fcsT^^ / dr [\p,(rf - ^ \p,{r)A (10) 



which is the first nonlinear correction, with 
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being the linear response counterion density orbital 
around a single macroion. In Eq. (p]), k can be identi- 
fied now as an inverse Debye screening length. The first 
and second terms on the right side of Eq. (^ account, re- 
spectively, for the counterion entropy and the macroion- 
counterion electrostatic interaction energy, while Eq. ^ 
corrects the latter term for nonlinear screening. 

The two- and three-body terms in Uejf involve effec- 
tive pair and triplet interactions. The pair interaction, 
w"'(r), is the sum of the linear response form § 
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i.e., the DLVO potential Q (aside from the excluded- 
volume adjustment to k), and the nonlinear correction 
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Finally, the triplet interaction is 



i;(3)(ri,r2,r3) - - fcsT^ / drpo(|ri-r|) 

X po(|r2-r|)po(|r3-r|). (14) 

Note that, with increasing salt concentration or decreas- 
ing macroion volume fraction {ndn^ — s- 0), the nonlinear 
corrections, Ai? and Ati(2)(r), and the effective triplet 
interaction, all diminish in relative magnitudes. Thus, 
nonlinear effects are strongest for deionized suspensions. 

Equations (0), and (|l|) are the main theoreti- 

cal results of the paper. Equation ([l3|), which is a hulk 
nonlinear correction to Ug (r) , is similar in structure to a 
wa/Z-mdwcerf pair interaction derived in Ref. jl^ [Eq. (13) 
therein]. Furthermore, Eq. (p^ is consistent with the 



triplet force derived in Ref. |13] [Eq. (33) therein] . The re- 
sponse theory, in neglecting short-range correlations be- 
tween microions, is formally equivalent to the mean-field 
PB theory. Its scope is therefore limited to long-range in- 
teractions between macroions. However, whereas numer- 
ical solutions of the nonlinear PB equation are practical 
only within symmetric Wigner-Seitz cells, characteristic 
of crystals (PB cell models) [||, the effective interac- 
tions derived here apply, in principle, to any bulk phase. 
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FIG. 1: Effective pair interaction between macroions of di- 
ameter a = 100 nm, valence Z = 700, and volume fraction 
ri = 0.01 at ionic strength Ca = 1 fJ,M. Solid curve: nonlin- 
ear screening. Dashed curve: linear screening (DLVO). Inset: 
cross-over from e"''' /r to e"""^ behavior at longer range. 



We can now quantify the practical significance of non- 
linear screening by evaluating the effective pair interac- 
tion and volume energy for selected system parameters. 
All results presented are for bulk, aqueous suspensions 
{Xb — 0.714 nm) and monovalent countcrions (z — 1). 
Analytic reduction of Eqs. (^0|) and (jl^) was checked 
by Monte Carlo integration. Figure 1 shows the effec- 
tive pair potential for macroions of diameter a — 100 
nm, valence Z = 700, and volume fraction r/ = 0.01, at 
ionic strength Cs = 1 /iM. The chosen charge is com- 
parable to the maximum renormalized charge (allow- 
ing for counterion condensation), Z* = Cu/Xb, where 
C — 0(10) At short and intermediate ranges, non- 
linear screening sharpens the counterion profile around 
a macroion, which, as seen in Fig. 1, weakens the pair 
repulsion. This result is consistent, if only qualitatively, 
with bulk measurements of pair interactions shorter 
in range than predicted by DLVO theory. At longer range 
(inset), nonlinear screening changes the asymptotic form 
to v^'^\r) ~ e"'^'', a qualitative departure from the more 
rapidly decaying DLVO form, Vq (r) ~ e'^"^ jr. As the 
macroion charge increases, higher-order nonlinear effects 
ultimately manifest themselves. In fact, for values of 
Z > Z* the effective pair potential can actually become 
attractive. Since pair attractions are mathematically ex- 
cluded within the PB theory [p^ , higher-order nonlinear 
terms must then also contribute in such extreme cases. 
Within physical parameter ranges (Z < Z*), however, 
?;'^^^(r) is always repulsive. Interestingly, the short-to- 
intermediate-range part of the nonlinear pair potential 
may still be reasonably fit by a DLVO-like potential. 

The cleanest test of the theory is comparison with 
ab initio simulation p^, which, like our response the- 
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FIG. 2: Total interaction potential energy for two macroions 
of diameter a = 106 nm and valence Z = 200, at zero ionic 
strength, in a cubic box of length 530 nm with periodic bound- 
ary conditions. Solid curve: nonlinear screening prediction. 
Dashed curve: linear screening (DLVO) prediction. Symbols: 
ab initio simulation data llTtl. 



ory, ignores short-range counterion correlations. Fig- 
ure 2 presents the comparison for the total potential en- 
ergy of interaction between a pair of macroions, of di- 
ameter a = 106 nm and valence Z = 200, in a cubic 
box of length 530 nm with periodic boundary conditions 
(taking into account image interactions) in the absence 
of salt. The agreement is essentially perfect, albeit for 
a case of relatively weak nonlinearity. Further simula- 
tions of more highly charged macroions are needed to 
more severely test the theory. A measure of the impact 
of nonlinear screening on thermodynamics is the rela- 
tive magnitude of Aw(^)(r) at the mean nearest-neighbor 
separation, r„„, where z;(^^(r) makes its dominant con- 
tribution to the potential energy. Figure 3 maps out, 
in the space of volume fraction and ionic strength, the 
boundary of the region within which | Aw^^^ (?"««) | exceeds 

(2) 

Vq {rnn) by at least 20% for the fee crystal structure: 
r„„/cr = 2-i/2(27r/377)i/3. With increasing Z and de- 
creasing Cs, the threshold 77 decreases. 

Thermodynamic properties of charged colloids are in- 
fluenced also by the volume energy, E, which contributes 
to the free energy, F, and pressure, P = —{dF/dV)T- 
The density dependence of E has been predicted - within 
linear screening - to drive an instability toward phase 
separation at low ionic strengths 0, The impact of 
nonlinearity is illustrated in Fig. 4, which plots the vol- 
ume energy contribution, P„ — -~{dE/dV)T, to the total 
equation of state for a system characterized hy a — 500 
nm, Z = 4000, and = 1 /xM. 

Differences between the linear and nonlinear predic- 
tions for P„ grow with increasing Z and decreasing c^, 
demonstrating the capacity of nonlinear screening to sig- 
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FIG. 3: Nonlinear deviations from DLVO theory for macroion 
diameter a = 100 nm and valences, from top to bottom, 
Z = 500, 600, 700. Systems with macroion volume fraction, 
rj, and ionic strength, Cs, above the curve deviate from the 
DLVO pair potential by at least 20% at the fcc-crystal nearest- 
neighbor distance. 
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FIG. 4: Volume energy contribution, Pv, to the total reduced 
pressure vs. volume fraction, rj, for macroion diameter a = 
500 nm and valence Z = 4000 at ionic strength Cs = 1 fj,M. 
Solid curve: nonlinear screening prediction. Dashed curve: 
linear screening (DLVO) prediction. 



nificantly affect phase stability of highly-charged, deion- 
ized suspensions. In fact, the parameter ranges in which 
instabilities toward phase separation have thus far been 
predicted |0, |8j lie well within the nonlinear regime. Fu- 
ture work will further explore influences of nonlinearity, 
including the role of attractive triplet interactions, on 
thermodynamics of charged colloids. 

In summary, applying response theory to the primitive 
model of charged colloids, we have consistently extended 
the DLVO linear-screening theory by including leading- 
order nonlinear contributions to effective interactions be- 



tween macroions. A major conclusion is that nonlinear 
corrections, aside from generating effective many-body 
interactions, can significantly modify the effective pair 
interaction and volume energy, and thus thermodynam- 
ics, of deionized suspensions of highly charged macroions. 
While mean-field pair interactions are repulsive at long 
range, nonlinear screening weakens shorter-range pair re- 
pulsion and generates effective three-body attractions. 
Nonlinear effects may thus facilitate any fluctuation- 
mediated short-range attractions p9| and help to explain 
anomalous phase behavior observed in experiments ^ , 
as well as in recent simulations | |20[ |. 

Discussions with C. N. Likos, H. Lowen, and 
K. S. Schmitz, and support of the Natural Sciences 
and Engineering Research Council of Canada (Grant 
No. 217232-1999), are gratefully acknowledged. 
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